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Recent developments in dynamically manipulating nanosystems have been drawing 
fresh and general interest in thermostatistics in a small scale. From the physics 
viewpoint, they seem to require deeper understandings of nonequilibrium 
thermostatistics. In this direction, some intriguing discussions have been made in the 
literature. Among others, an attempt of Jarzynski [1] to connect nonequilibrium works 
with equilibrium free energies has been attracting much attention both theoretically [2] 
and experimentally [3]. It leads to necessity of reexamining the basic thermodynamic 
quantities in the nonequilibrium regime. 

In this paper, we study the concept of work in nonequilibrium quantum 
thermostatistics. We find that there naturally emerges the gauge structure in the theory, 
the field of which is referred to as the work gauge field. In particular, we consider a 
cyclic process in the space of the parameters, which describe the external driving of the 
system. Then, the work turns out to be given by the flux of the work gauge field. We 
examine this general result for an example of a spin- 1/2 interacting with an external 
magnetic field. We also make comments on a possibility of observing such a geometric 
effect, for example, in an NMR experiment and its relevance to the problem of 
cooling/heating of a small system. 

Let us start our discussion with the quantum mechanical formulation of the first law 
of thermodynamics. Consider a driven quantum system in contact with an environment. 
Its state is represented by a density matrix p, which is a Hermitian, traceclass, positive 

2 



semidefinite operator. The internal energy is given by the expectation value of the 
system Hamiltonian H with respect to p: U = Tr (p H) . Its infinitesimal change is 
given by 



dU = Tr (dp H) + Tr (p dH) . (1) 

Then, the first law of thermodynamics, d' Q- dU + d'W , is realized if the following 
identifications are made for the work and the quantity of heat: 

d'W =-Tr(p dH), (2) 

d'Q = Tr(dpH). (3) 

The external driving is described by the dependency of the Hamiltonian on a set of 
the parameters 

H = H (A 1 , A 2 ,--, A"). (4) 

Note that, in contrast to Ref. [1], here we are considering the multiple parameters. 
Accordingly, the work in Eq. (2) is written as 



d' W = Y, a^dV, (5) 



where 



a„=-Tr(pd M #) 



(6) 



with the notation, d ^ = d I d^. As we shall see below, the quantity in Eq. (6) is a 
gauge field potential, which is referred to here as the work gauge field. It may be 
analogous to the gauge field in the discussion of the Berry phase [4] in the sense that it 
is induced by the external parameters. 

Now, for a process along a closed curve C in the parameter space, the work is 
written as 

W c = £§a^dV=l£ \\f^ v dVAdr. (7) 

M=l C ^ fi,v=l s 

where / = d ^ a v - d v and C-dS(i.e., C is the boundary of the surface 5 in the 
parameter space). 

Clearly, the "field strength", / , vanishes (and thus W c =0) if the density matrix 
is a function only of the Hamiltonian. Therefore, for example, the equilibrium state 
represented by the canonical density matrix 

with the partition function, Z(f5) = Tr e~^ H , does not yield the nonvanishing field 
strength. In fact, for p eq , the work gauge field becomes pure gauge, a^=-d , 
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where F = -/5~ 1 lnZ(/3) is the free energy. This is nothing but a familiar relation that 
the work done is given by the free energy difference in an isothermal process. 

The above discussion tells us that the associated gauge transformation of the density 
matrix is 

p^p + p, (9) 

where p is a traceless matrix depending only on the Hamiltonian and does not violate 
the positive semidefiniteness of the transformed density matrix. Under Eq. (9), the work 
gauge field changes as follows: 

a^a^+d^A, (10) 

where d ^ A = -Tr [p(H) d ^ H] (such A exists since p is a function only of the 
Hamiltonian). is, therefore, an Abelian gauge field. The work in Eq. (7) is a gauge- 
invariant quantity. 

Thus, the work gauge field can give rise to a nonvanishing field strength if the state 
is out of equilibrium, in general. To realize such a nontrivial work gauge field, we 
employ a quantum operation [5] on the equilibrium canonical density matrix in Eq. (8). 
That is, p eq is linearly mapped to a nonequilibrium state p as follows: 

p eq ^p=$(p eq ). (11) 



5 



In particular, we represent the map by using the positive operator-valued measure 
(POVM) [5] 

P = *(P eq ) = L V kPeq V k l (12) 

k 

It satisfies the trace-preserving condition, V^ 1 V k = I, where / stands for the 
identity operator and V k 's may depend on the set of the parameters {A,} f=1 2 ... n . To 
obtain a nonequilibrium state, clearly V k 's cannot be functions only of the Hamiltonian. 
Unlike in the equilibrium state p eq , f5 does not have the meaning as the inverse 
temperature in the nonequilibrium state p, in general. Substituting Eq. (12) into Eq. (6), 
we obtain the nonvanishing work, W c . 

We wish to notice two point, here. One is that there are a variety of nonequilibrium 
states, and the present construction is, though quite general, just one example. However, 
the above one is useful since it is in reference to the equilibrium state. The other is that, 
due to the trace-preserving condition, a monopole-type singularity does not appear in 
the configuration of the work gauge field. 

Finally, let us consider an example of a single spin- 1/2 interacting with a magnetic 
field B = (B x , B y , B z ), which plays a role of the external parameters. The system 
Hamiltonian reads 
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H = k<jB. (13) 

Here a ' s are the Pauli matrices and K is a constant involving the magnetic moment. 
The Planck constant as a multiplicative factor is set equal to unity for the sake of 
simplicity. The equilibrium canonical density matrix is given by 



Peq 



/-^-^tanh(jBjrfl) 
B 



(14) 



where B = | B | . 

To map the equilibrium state in Eq. (14) to a nonequilibrium one, we consider the 
following most general POVM: 

V k =c k I + X k -a, (15) 

where c k 's and X/s are complex numbers and complex vectors generically 
dependent on the magnetic field, respectively. The trace-preserving condition leads to 

W|c,| 2 +X;-X,) = l, (16) 

k 

j;( c ;x t+ c t x;+ix;xxj=o. a?) 

k 

Using the nonequilibrium density matrix p = V k p eq V/ with Eqs. (14) and (15), we 
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find that the work gauge field is given by 



HB) = -K 1 £(c k Xl + c;X k + iX k xXl) 



.tanh(jB kB) 



+ k — -> | Icl B + /c t X t xB + /c t BxX 



+ (B-X t )X;+(B-X;)X t -(X;-X t )B]. (18) 
This field is drastically simplified if the following choice is made: 

X t = (0,0,* t ), (19) 

where g k 's , as well as c k 's, are real constants. In this case, the condition in Eq. (16) 
and (17) become 

EK 2 + * t 2 ) = l, (20) 

k 

l,c kgk = 0, (21) 

k 

respectively. Now, performing a straightforward calculation, we obtain the work gauge 
field of the following form: 

a(B)= f(B)(aB x ,ccB y ,B z ), (22) 



where 



f(B)=K 



tanh(ft kB) 
B 



(23) 



(24) 



k 



Taking the rotation of this field, b(B) = V B x a(B) with the notation V B = d/dB, we 
have the following field strength: 



where is the unit vector in the direction of rotation around the 5 z -axis, that is, 



For a process along a closed curve C in the B-space, the work is given in terms of 
the flux of b traversing the surface S surrounded by C. In other words, given the 
work gauge field strength in Eq. (23), the work is determined sorely by the geometric 
configuration of C . 

In a particular case when all c k 's vanish, a = -1 and the quantum operation 
becomes isoentropic: 4>(p eq ) = o z p eq o z . 

In conclusion, we have studied the concept of work in quantum thermostatistics of a 




(25) 




(26) 
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system driven by an external force. We have found that the gauge theoretical structure 
naturally emerges in a nonequilibrium process, which is generated by using the positive 
operator-valued measure (POVM). We have described the thermodynamic work as the 
flux of the gauge field termed the work gauge field and have considered it for a cyclic 
process in the space of the external-force parameters. The example of a spin- 1/2 
interacting with an external magnetic field has been explicitly analyzed. It is our hope 
that such an effect can be observed, for example, in an NMR experiment. 

Since the external parameters and the POVM operators are experimentally 
controllable, one can, for example, reverse the direction of the cyclic process C. Under 
such a reversal, W c changes its sign. Therefore, the present study can be applied to the 
problem of cooling/heating of a small system, which is of physical relevance to 
nanoscience. 
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